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By using the concept of negativity, we investigate entanglement in (1/2,1) mixed-spin Heisenberg 
systems. We obtain the analytical results of entanglement in small isotropic Heisenberg clusters 
with only nearest-neighbor (NN) interactions up to four spins and in the four-spin Heisenberg model 
with both NN and next-nearest-neighbor (NNN) interactions. For more spins, we numerically study 
effects of temperature, magnetic fields, and NNN interactions on entanglement. We study in detail 
the threshold value of the temperature, after which the negativity vanishes. 
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I. INTRODUCTION 

Entanglement, an essential feature of the quantum me- 
chanics, has been introduced in many fields of physics. In 
the field of quantum information, the entanglement has 
played a key role. The study of entanglement properties 
in many-body systems have attracted much attention 0]- 
psf. The Heisenberg chains, widely studied in the con- 
densed matter field, display rich entanglement features 
and have many useful applications such as in the quan- 
tum state transfer [2^ . 

Most of the systems considered in previous studies are 
spin-half systems as there exists a good measure of entan- 
glement of two spin- halves, the concurrence |30| . which 
is applicable to an arbitrary state of two spin halves. 
On the other hand, the entanglements in mixed-spin or 
higher spin systems are not well-studied due to the lack 
of good operational entanglement measures. There are 
several initial studies along this direction [sil Is^ l33l |. 
however these works are restricted to the case of two 
particles. 

For the case of higher spins, a non-entangled state has 
necessarily a positive partial transpose according to the 
Peres-Horodecki criterion |Q. In the case of two spin 
halves, and the case of (1/2,1) mixed spins, fortunately, 
a positive partial transpose is also sufficient. Thus, the 
sufficient and necessary condition for entangled state in 
(1/2,1) mixed spin systems is that it has a negative par- 
tial transpose. This allows us to investigate entanglement 
features of the mixed spin system. 

The Peres-Horodecki criterion give a qualitative way 
for judging whether the state is entangled or not. The 
quantitative version of the criterion was developed by 
Vidal and Werner j35 |. They presented a measure of 
entanglement called negativity that can be computed ef- 
ficiently, and the negativity does not increase under local 
manipulations of the system. The negativity of a state p 
is defined as 



Mi 



(1) 



where /Zi is the negative eigenvalue of p'^^ , and Ti denotes 
the partial transpose with respect to the first system. 



The negativity Af is related to the trace norm of p-^^ via 

^■W = ^^. (2) 

where the trace norm of is equal to the sum of the ab- 
solute values of the eigenvalues of p^^ . In this paper, we 
will use the concept of negativity to study entanglement 
in (1/2,1) mixed-spin systems. 

As shown in most previous works, models with the 
NN exchange interactions are considered and it is not 
easy to have pairwise entanglement between the NNN 
spins 0, 0|. It is true that there exist some quasi- 
one-dimension compounds offering us systems with NNN 
interactions. Bose and Chattopadhyay (3^] and Gu 
et al. 37] have investigated entanglement in spin-half 
Heisenberg chain with NNN interactions. In our paper 
here, we study entanglement properties not only in the 
(1/2,1) mixed-spin systems only with NN interactions, 
but also in the system with NNN interactions. 

Entanglement in a system with a few spins displays 
general features of entanglement with more spins. For 
instance, in the anisotropic Heisenberg model with a 
large number of qubits, the pairwise entanglement shows 
a maximum at the isotropic point |3^ . This feature 
was already shown in a small system with four or five 
qubits [ill- So, the study of small systems is meaningful 
in the study of entanglement as they may reflect general 
features of larger or macroscopic systems. Also, due to 
the limitation of our computation capability, we only con- 
centrate on small systems such as 4, 5 and 6-spin models. 

The paper is organized as follows. In Sec. II, we study 
the systems with only NN interactions. The analytical 
results of negativity for the cases of two and three spins 
are given. The relation between entanglement and the 
macroscopic thermodynamical function, the internal en- 
ergy is revealed. Also we numerically compute the neg- 
ativity in more general mixed-spin models up to eight 
spins, and consider the effects of magnetic fields in this 
section. In Sec. HI, the system with NNN interaction 
is discussed. For the four-spin case, we analytically cal- 
culate the eigenenergy of the system from which we get 
the analytical results of the negativity of the NN spins. 
We numerically study negativities versus NNN exchang- 
ing coupling, and the case of finite temperature is also 



considered. For larger system up to eight-spin system, 
we get some numerical results. The conclusion is given 
in Sec. IV. 



II. ENTANGLEMENT IN HEISENBERG CHAIN 
ONLY WITH NEAREST-NEIGHBOR 
INTERACTION 

A. Analytical results of Entanglement in 
Heisenberg models 

We study entanglement of states of the system at 
thermal equilibrium described by the density operator 
p{T) = exp{-PH)/Z, where /3 = l/ZcsT, fcs is the Boltz- 
mann's constant, which is assumed to be one throughout 
the paper , and Z = Tr{exp(— is the partition func- 
tion. The entanglement in the thermal state is called 
thermal entanglement. 

We consider two kinds of spins, spin ^ and 1, alternat- 
ing on a ring with antiferromagnetic exchange coupling. 
The Hamiltonian is given by 



N/2 

Ho = ■ S, + S, ■ s,+i), (N e even) (3) 

Hq = (Sj ■ Sj + Sj ■ Sj+i) + S (jv + i) ■ Si, 

i=l 

{N e odd), (4) 

where and Si are spin-1/2 and spin-1 operators, re- 
spectively. The exchange interactions exist only between 
nearest neighbors, and they are of the same strength 
which are set to one. We adopt the periodic boundary 
condition. In Fig. 1, we give the schematic representa- 
tion of the above Hamiltonian. Next, we first consider 
the models with two and three spins, and aim at getting 
analytical results of entanglement. 



1. Two-spin case 

For the two-spin case, the Hamiltonian Q reduces to 
Hq = Si • Si. To have a matrix representation of the 
Hamiltonian, wc choose the following basis 

{l-^, -1), 1^,0), 1-1, 0),|i,l>}, (5) 

where |m, M) is the eigenstate of Sz and Sz with the cor- 
responding eigenvalues given by m and M, respectively. 

In the above basis, the Hamiltonian can be written as 
a block-diagonal form with the dimension of each block 
being at most 2x2. Thus, the density matrix pi2 for the 
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FIG. 1: Schematic representation of the mixed spin Hamilto- 
nian for the number of spins from 2 to 5. 



thermal state is obtained as 
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with the partition function and the matrix elements given 
by 

4/3 



z 


= 2e^-\ 


h 4e 


a\ 


= 06 = 


1 

e 2 
1 


a2 


= 05 = 


6' 


as 


= 04 = 


K 

V2 


bi 


= b2 = 


3 



/z, 



/z, 
/z, 



V2{ai — 02). 



(7) 
(8) 
(9) 

(10) 



(11) 



After the partial transpose with respect to the first spin- 
half subsystem, we can get 
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(12) 



which is still of the block-diagonal form, and computation 
of its eigenvalues is straightforward. There are only two 
eigenvalues which are possibly negative. The negativity 
is thus given by 



2 
1 
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2 



J^ipu) =^ max [0, (fli - 02)2 -t- 4bl - m - 02)] 

[0,^(05-06)2 + 452-05-06]. (13) 
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Substituting Eqs. lf7|l- Hll() leads to the analytical result 
of negativity 



N{p) = max [0, ^J (ai - az)^ + Abl 
—2 max[0, 02 — 2ai] 



= — max 
3 



0,- 



4e- 



;/3 + 2e-5/3 



(14) 



where the second equality follows from Eq. Hll() . 

We can see that the negativity is a function of the 
single parameter /3. In the limit of T ^ 0, the nega- 
tivity becomes 1/3 and the ground state is entangled. 
From Eq. H14|l . it is direct to check that the negativity 
is a monotonically decreasing function when temperature 
increases. After a certain threshold value of the tempera- 
ture, the entanglement disappears. This threshold value 
Tth can be obtained as 



Tth = 3/(41n2) w 1.0820. 



(15) 



For a ring of spin-half particles interacting via the Hiscn- 
berg Hamiltonian, it was shown that the pairwise thermal 
entanglement is determined by the internal energy [3^ . 
It is natural to ask if similar relations exist in the present 
mixed-spin system. The internal energy can be obtained 
from the partition function as 



Z d(3 



(16) 



Substituting Eq. Q into the above equation leads to 



u = 



e'^ + 2e- 



(17) 



From Eqs. (|14|l and H17|) . we obtain a quantitative rela- 
tion between the negativity and the internal energy 



Af{p) = - max[0, - 



4- 2C/]. 



(18) 



The above equation builds a connection between the 
microscopic entanglement and the macroscopic thermo- 
dynamical function, the internal energy. The internal 
energy completely determine the thermal entanglement. 
From the equation, we can also read that the thermal 
state becomes entangled if and only if the internal en- 
ergy U < -1/2. Since U = {H) = (si • Si), we have 



AA(p) - imax[0,-l-2(si -Si)], (19) 

which is consistent with the result obtained in Ref. jsj 
by the group-theoretical technique. 



2. Three-spin case 

We now consider the three-spin case and the schematic 
representation of the corresponding Hamiltonian is given 



by Fig. 1. In this situation, there are two types of pair- 
wise entanglement, the entanglement between spin 1/2 
and spin 1 and the entanglement between two spin halves. 

The eigenvalue problem can be solved analytically, and 
after tracing out the third spin-half system the reduced 
density matrix pi2 is still of the same form as in Eq. © 
with matrix elements given by 



fli 



02 



as 




Z = 



(20) 



Substituting the above equations to Eq. (|13|l leads to the 
negativity 



M{pi2) = max 



1/5 



/Z. (21) 



It is evident that the negativity becomes 1/3 in the limit 
of T ^ 0. From the expression of the negativity, the 
threshold value of temperature after which the entangle- 
ment vanishes can be estimated as 



Tth ~ l/ln3.2719« 0.7609. 



(22) 



To examine the entanglement between two spin halves, 
we trace out the spin 1 system and get the reduced den- 
sity matrix pi3 as follows 
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(23) 



with the matrix elements given by 
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After taking the partial transpose, we can get pj^ 
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Then, the negativity is readily obtained as 
Af{pi3) =max[0, l&l -oi]. 



(25) 



(26) 
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It is straightforward to check that the negativity is always 
zero. Or, from another way, aU the eigenvalues of the 
matrix pf^ are obtained as 

A2 = (^le'if' + 2e'^^ + lei^^/Z, 

A3 - A4 = a2. (27) 

Obviously the negativity vanishes here, in other words 
there is no entanglement between the two spin halves. 

The ground-state negativity A/12 = A/'(pi2)=l/3 and 
A/13 — A/'(pi3)=0, here A/12 denotes the negativity be- 
tween the 1th and 2th spin on the chain and A/13 denotes 
the one between the 1th and 3th spin. The equation above 
can be obtained from the non-degenerate ground state 
given by: 

It is interesting to see that the ground-state entanglement 
between the spin half and spin 1 in the three-spin case is 
the same as that in the two-spin case. 

Due to the SU(2) symmetry in our system, there are 
following relations between correlation functions and neg- 
ativities 

A/'i2=-^-^(si-Si), 

AA23=-i-^(Si.S2), 

AA3i=-i-(s2-si), (29) 

where we have removed the max function in the negativ- 
ity, implying that the negative value of Af indicates no 
entanglement. Then, we have the relation between the 
internal energy and the negativities 

5 3 

[/ = --- AA13 - -(A/'i2 + A/'23) 

= -^-AAi3-3AAi2. (30) 

The second equality follows from the exchange symmetry, 
namely, the Hamiltonian is invariant when exchanging 
two spin halves. So, for the three-spin case, the internal 
energy is related to two negativities. 

To apply the above result, we consider the the ground- 
state properties (T = 0). The Hamiltonian can be rewrit- 
ten as 

i^ = ^[(si + Si+S2)2-s?-S?-s2]. (31) 



Then, by the angular momentum coupling theory, the 
ground-state energy is obtained as Eq = —7/4. Sub- 
stituting the ground-state energy and A/12 = 1/3 to 
Eq. H3U|I . we obtain A/13 = ~l/2, indicating that there 
exists no entanglement between two spin halves. Next, 
we consider more general situations, i.e., the case of even 
N sites. 

3. The case of even N spins 

Except for the SU(2) symmetry in the system, there 
exists exchange symmetry for the case of even spins. For 
instance, in the four-spin model, the Hamiltonian is in- 
variant when exchanging two spin halves or two spin ones. 
Thus, for even-spin model, the entanglements between 
the two nearest-neighbor spins and the correlation func- 
tions {si ■ Si) are independent on index i. Therefore, the 
internal energy per spin is equal to the correlation func- 
tion (sj • Si) 

u^U/N ^ {si- Si). (32) 
From Eqs. and we have 

^2 = -^-^". (33) 

This equation indicates that for the case of even spins 
the entanglement between two nearest neighbors is solely 
determined by the internal energy per spin. And for the 
case of zero temperature, the entanglement is determined 
by the ground-state energy. The less the energy, the more 
the entanglement. 

We now apply the above result to the study of ground- 
state entanglement. We rewrite the four-spin Hamilto- 
nian as follows 

H = i[(si + Si+S2 + S2)2-(si+S2)2-(Si+S2)2]. (34) 

Then, by the angular momentum coupling theory, we ob- 
tain the ground-state energy per site eq = i?o/4 = —3/4. 
Thus, from Eq. we have Afi2 = 1/6. For > 5, it is 
hard to get analytical results. We next numerically cal- 
culate the entanglement for the case of more spins, and 
also consider the effects of magnetic fields. 

B. Numerical results 

Having obtained analytical results of entanglement in 
the Heisenberg model with a few spins, we now nu- 
merically examine the entanglement behaviors in more 
general Hamiltonian including more spins and magnetic 
fields. 

1. Entanglement versus temperature 

We consider the entanglement versus temperature for 
different number of spins N, and the numerical results 
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2. Effects of magnetic fields 




FIG. 3: Negativity versus temperature for different odd num- 
ber of spins in the Heisenberg model. 



are plotted in Fig. 2 and Fig. 3. It is clear to see that the 
ground state exhibits maximal entanglement, and with 
the increase of temperature, the entanglement monoton- 
ically decreases until it reaches zero. The decrease of en- 
tanglement is due to the mixture of less entangled excited 
states when increasing the temperature. The existence of 
the threshold temperature is also evident. For the case of 
even number of sites (Fig. 2), the threshold temperature 
decreases with the increase of A'^. In contrast to this be- 
havior, for the case of odd number of spins, as seen from 
Fig. 3, the threshold temperature increases with the in- 
crease of N. 



We now examine the effect of magnetic fields on entan- 
glement. The Heisenberg Hamiltonian with a magnetic 
field along z direction is given by 



Hi — Hn 



N/2 

{stz+ S^,^),{N e even). 



(35) 



i=l 



For N = 2, the analytical result of negativity can be 
obtained (see Appendix A). 

In Fig. 4, we plot the negativity versus the magnetic 
field at a low temperature. For the two-spin case, with 
the increase of the magnetic field, the negativity rapidly 
reach a platform, and then after a certain magnetic field 
Bth, it jumps down to zero, indicating no entanglement. 
When B — 0, the ground state is two-fold degenerate 
with TV = 1/3. For a small B > 0, the system is no 
longer degenerate, and the ground state is given by 

IV^o) = ^|i-l)-f |~i0), (36) 

which is of the Schmidt form. For a state written as its 
Schmidt form 



1^) =^c„|e„)|/„). 



(37) 



the negativity is obtained as |3£ 



- 1 



/2. 



(38) 



Here, C„ are the Schmidt coefficients, and {|e„)} and 
{!/„)} are bases for subsystems 1 and 2, respectively. 
Applying this formula to the ground state, we imme- 
diately have Af = When B = Bth = 3/2, the 
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ground state is two-fold degenerate, and when B > i3th 
the ground state becomes non-degenerate and the corre- 
sponding wave function is given by 

IV'o) = 1-^,-1), (39) 

which is obviously of no entanglement. Then, the ground- 
state negativity forms a platform when < B < Bth- 
The jump of negativity a,t B = Bth is due to the level 
crossing. For N > 2, the effects of magnetic fields on 
entanglement can be also explained by level crossing. For 
instance, for = 4, there are two level crossing, and the 
entanglement displays two jumps. 



III. EFFECTS OF NEXT-NEAREST-NEIGHBOR 
INTERACTIONS ON ENTANGLEMENT 

We have studied the effects of finite temperature and 
magnetic fields on entanglement, and now consider the 
model containing two kinds of spins, spin ^ and 1, alter- 
nating on a ring with antiferromagnetic exchange cou- 
pling between both the NN spins and the NNN spins. 
The Hamiltonian can be expressed as 

N/2 N/2 

H = Ji^^ (si-Si + Si-Sj+i) + J2 ^ (sj-Si+i + Si-Si+i), 

i=l i=l 

(40) 

where the and Si are spin- 1/2 and spin-1 operators in 
the ith cell. Ji characterizes the NN exchange coupling 
and J2 the NNN coupling. We consider the antiferro- 
magnetic interaction by taking Ji, J2 > 0. N is the total 
number of spins and here we choose it be even. Also we 
adopt the periodic boundary condition. 



A. Four-spin model 

1. Eigenenergy and ground-state entanglement 

The model with four spins is the simplest model with 
NNN interactions. We first solve the eigenvalue problem 
of this model. The key step is to write the four-spin 
Hamiltonian in the following form, 

H = i{JiS2 + (2J2- Jl)[(si+S2)2 + (Sl+S2)2] 

- 2Msl + sl + Sl + Sl)}, (41) 

where S = Si -|- S2 + Si -I-S2 denotes the total spin. From 
the above form, by angular momentum coupling theory, 



one can readily obtain all eigenvalues of the system as 
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where parameter n — 0, 1, 2 in expressions Ei and E2 and 
n = 1,2 in expressions £^3 and i?4, respectively. Then, 
from Eq. 142|) . we may find the ground-state energy as 

r -3Ji + I J2 when < J2 < Ji/4, 
£^GS = < -2Ji - I J2 when Ji/4 < J2 < ^i/2, 
[ -^J2 when J2 > Ji/2, 

(43) 

Clearly, there are two level-crossing points, which will 
greatly affect behaviors of ground-state entanglement. 

From Eq. (|4HI . it is obvious that in addition to the 
the SU(2) symmetry, there also exists an exchange sym- 
metry, namely, exchanging two spin halves or two spin 
ones yields invariant Hamiltonian. Then, the correlator 
(si • Si) between any NN spins are the same. Therefore, 
we can get the correlator (s.^ • S^) from the ground-state 
energy via the Hellmann-Feynman theorem. 

When < J2 < Ji/4, after applying the Hellmann- 
Feynman theorem to the ground-state energy, we obtain 

Due to the SU(2) symmetry in our system, we have the 
following relation between the negativity and the corre- 
lator (si • Si)[3l|, 

AAi/2,i=max{0,-i-^(si.Si)}, (45) 

where we use A/1/2,1 to denote the negativity between the 
NN spins. Thus we obtain 

Ml/2.1 = 1/6. (46) 

In other regions, we find that A/1/2,1 = for J2 > Ji/4. 

From the above analytical results, we can find that 
the negativity between the NN spins is not a continuous 
function of parameter J2. It jumps from the value ^ down 
to zero at J2 = 0.25. Hence, we can see that the NNN 
interaction may deteriorate the entanglement of the NN 
spins. The critical point of J2 — 0.25 may be considered 
to be a threshold value, after which the NN entanglement 
vanishes. 

Having studied ground-state entanglement, now we 
make a short discussion of entanglement of excited states. 
We consider the first excited state. In the region Ji/4 < 
J2 < 3Ji/8, the first excited energy is given by 

El = -3Ji -f ^ J2, (47) 
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FIG. 5: Negativity versus the exchange interaction J2 at the 
temperature T = 0.008 in the four-spin system. The ground- 
state and the first two excited-state energy levels versus J2 
are inserted. 



from that we can get (si • Si) = —3/4 and -/Vi/2,1 — 1/6. 
While in the rest region the negativity of the first excited 
state is zero. 



2. Thermal entanglement 

Next, we consider the thermal entanglement. From 
Eq. H42|) . the partition function can be obtained as 



e 2 



-I- 5e 



3e 



3/3 Ji 



(48) 



The correlator (si - Si) at finite temperature can be com- 
puted from the partition function via the following rela- 
tion 



(si - Si) =-. 



1 dZ 



4:(3Z dJi ' 
Substituting lO to Eq. (gHJ yields 



(49) 



(si • Si) = - — [e-i'^J^-Ue-^l'J' + he^-'^ -f Qe"^-'^) 

+ ei^^^i^be-^-^^ +?,e'^'^^ +2e^'^-^% (50) 

After substituting the above equation into (|45|) . we may 
get analytical expression of the negativity A/1/2.1 ^-t finite 
temperatures. The negativity is a function of Ji, J2 and 
T. 

Low-temperature case: We now make numerical cal- 
culations of entanglement and first consider the low- 
temperature case. We take Ji = 1 in all the following 
plots. In our system there exist three kinds of negativ- 
ity, the negativity Ai/2,1 between spin-1 and spin-1/2, 



A/1/2,1/2 between two spin-1/2, and A/1,1 between two 
spin- 1 . 

In Fig. 5, we plot the negativity versus J2 in four-spin 
system at a low temperature of T = 0.008. It is clear 
to see that A/1/2,1 keeps a value about 1/6 when J2 in- 
creases from zero until it reaches the critical point, at 
which the A/1/2,1 displays a jump to zero. This behavior 
of entanglement is consistent with that at zero tempera- 
ture from the analytical results. It is natural to see that 
increase of NNN exchange interaction will suppress the 
entanglement of NN spins, and at last completely erase 
the entanglement. 

In comparison with A/1/2,1, the negativities A/1,1 and 
•/Vi/2,1/2 behave distinctly. We see that near the point 
of J2 = 1/4, A/1,1 increases quickly to a steady value 
about 1/3, and when J2 reaches the value about 1/2, 
A/1,1 jumps another steady value 1. These two jumps re- 
sult from the two level crossing as seen clearly from the 
figure inserted. The second level crossing also leads to a 
small dip in the curve of A/1,1. The negativity A/1/2,1/2 
displays a jump to a steady value of 1/2 near J2 — 1/2. 
It is evident that the entanglement between NNN spins 
is enhanced by increasing NNN interactions. The com- 
petition between NN and NNN interactions leads to rich 
behaviors of quantum entanglement. Another observa- 
tion is that there is a range of J2, at which negativities 
A/1/2,1 and A/1/2,1/2 are zero, and only A/1,1 is not zero. 
This indicates that the NNN interaction must be strong 
enough to build up the entanglement of two spin halves. 

Entanglement versus J2 and T: As temperature in- 
creases the entanglement will decrease due to the mixing 
of less entangled excited states to the thermal state. It is 
obvious that there exists a threshold temperature after 
which the negativity is zero. In the frustrated system, 
there exists the parameter J2, and with its increase, the 
negativity A/1/2,1 will decrease to zero, while A/1,1 and 
A/1/2.1/2 increase from zero to their maxima. So it is clear 
that there also exists a threshold J2th corresponding to 
the boundary between zero and nonzero negativities. 

In Fig. 6, we show the negativity A/1/2,1 versus the 
temperature and J2. When the temperature approaches 
zero, the A/1/2,1 reaches its maximum, and then with the 
temperature increasing, A/1/2,1 decreases to zero. On the 
J2 — T plane, there is a curve along which the negativity 
just turns to be zero. It is possible to consider that the 
curve describes the threshold J2th versus the tempera- 
ture. Obviously, J2th does not behave as a monotonous 
function of the temperature, and it displays a peak at 
about T — 0.178, This behavior is in big contrast with 
the case of non- mixed qubit systems [s^ • When the tem- 
perature rises, the weight of excited states will increase 
and it may strongly affect the negativity. This behavior 
of J2th results from both the mixture of excited states 
to the thermal state and the intrinsic properties of the 
mixed-spin system. In addition, after crossing the tem- 
perature about T = 1.082, A/1/2,1 will vanish, irrespective 
of the value of J2- 

From another point of view, we can read the threshold 
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FIG. 6: Negativity A/1/2,1 versus both the exchange interac- 
tion coefficient J2 and the temperature in the four-spin sys- 
tem. 

temperature Tth versus different J2 from the curve in the 
J2 — r plane. When J2 increases, Tth decreases, and 
when J2 crosses about 0.3758, A/1/2,1 will disappear at 
any temperature. 

Next, we consider the entanglement between NNN 
spins. In Fig. 7, we plot the negativity A/1/2,1/2 as a func- 
tion of the temperature and J2. We can see that, before 
J2 reaches the value about J2 = 0.5, A/1/2.1/2 keeps being 
zero at any temperature. And in the region J2 > 0.5, the 
A/1/2,1/2 c'An be enhanced by the increasing NNN interac- 
tion. This is a result from the competition of two kinds of 
exchange interactions. The thermal fluctuation all along 
suppresses the entanglement. So, from the curve lying 
on the J2 —T plane which corresponds to the boundary 
of the nonzero and zero values of A/1/2,1/2, we may find 
that the higher the temperature is, the larger the thresh- 
old J2th will be. From another point of view, the Tth 
increases as J2 increases. 

In Fig. 8, we plot the negativity A/1.1 versus T and J2. 
In the region of J2 < 0.25, A/1.1 is zero at any temper- 
ature. When J2 > 0.25, the increasing NNN exchange 
interaction J2 enhances the negativity and exhibits two 
particular flat roofs. With the temperature rises, A/1,1 is 
suppressed to zero. Also we can consider the threshold 
Tth and J2th from the critical curve on the J2 — T plane, 
and the J2th also behaves as an increasing function of the 
temperature. We can see the nonzero region of A/1,1 is 
much larger than A/iy2, 1/2- But here it should be pointed 
out that, because A/1,1 > only gives a sufficient condi- 
tion for entangled state, we can not definitely say that 
the state in the area of zero negativity is not entangled. 



B. Numerical results of negativity for more spins 

In this section, we present numerical results of negativ- 
ity for more spins, and first consider the low-temperature 
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FIG. 7: Negativity A/1/2,1/2 versus both the exchange inter- 
action coefficient J2 and the temperature in four-spin system. 
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FIG. 8: Negativity A/1,1 versus both the exchange interaction 
coefficient J2 and the temperature in four-spin system. 



case. 



1. Low-temperature case 

In Fig. 9, we plot the negativity versus J2 for the case 
of six spins at a lower temperature. The negativity A/1/2,1 
behaves as a decreasing function of J2. It decreases to 
zero at about J2 = 0.27, which is the special point corre- 
sponding to the energy levc;l crossing. On the contrary, 
around the point J2 = 0.27, A/1,1 jumps up to a nonzero 
value, and then increases gradually until approaching the 
limit about A/i.i = 0.33. This behavior is quite different 
from that in the four-spin model. 

We also see that AAi/2,1/2 is zero all the time. It can 
be understood as follows. In the six-site system there are 
three spin halves with the NNN interaction. Even for a 
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FIG. 9: Negativity versus the exchange interaction J2 at the 
temperature T = 0.02 in six-spin system. The ground and 
the first two excited energy levels versus the J2 are inserted 
in the plot. 
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FIG. 10: Negativity versus the exchange interaction J2 at the 
temperature T = 0.02 in eight-spin system, The ground and 
the first three excited energy levels versus the J2 are inserted. 



pure homogeneous three-qubit system, there is no entan- 
glement between two spins, irrespective of the strength 
of the exchange interactions 5]. Now, in addition to the 
interaction among three spin halves, there are also inter- 
action between spin halves and spin ones. So, it is rea- 
sonable that the entanglement between two spin halves 
is zero. 

The negativity versus J2 in the eight-site case is shown 
in Fig. 10. After the first sharp jump to a value (not 
zero) at about J2 = 0.25, A/1/2,1 goes down to zero grad- 
ually. At approximately J2 — 0.55, the negativity is zero. 
On the contrary, the negativity A/1,1 jumps up at about 
J2 = 0.25, and then goes up gradually and almost lin- 
early until J2 reaches about 0.55. Then there happens 
a sharp decrease to nearly zero, and after that it begins 




FIG. 11: Negativity A/1/2,1 versus both the exchange interac- 
tion J2 and the temperature in the six-spin system. 




FIG. 12: Negativity 7Vi,i versus both the exchange interaction 
J2 and the temperature in six-spin system. 



to increase gradually. The negativity A/1/2,1/2 keeps zero 
until J2 reaches about 0.67, and then it goes up until 
reaches a steady value. 

In Fig. 10, we find that the three kinds of negativity 
exhibit different properties. From figure inserted, i.e., the 
energy levels of the eight-spin system, we can see that in 
the region from about J2 — 0.55 to J2 = 0.67 the first 
excited energy is quite close to the ground energy. It is 
known that the energy level crossing can greatly affect 
the entanglement. Here, the two close energy levels also 
play an important role in the behavior of entanglement. 
The approximate degenerate energy levels may remark- 
ably change the probability distribution even at a very 
low temperature, thereby affect the negativity. 
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2. Entanglement versus J2 and T 

Now, we present the entanglement versus J2 and T in 
the six-spin model. The NN negativity A'1/2,1 is shown 
in Fig. 11. On the J2 — T plane, similar to the four-spin 
case, the J2th does not behave as a monotonous function 
of T and it reaches its maximum at about T — 0.27. From 
the figure, we may find that the entanglement only exists 
in the region approximately T < 0.925 and J2 < 0.418. 
The strong NNN interaction and thermal fluctuation will 
suppress the NN entanglement to zero. 

We do not plot the negativity A/1/2.1/2 as it is zero all 
along for any J2 and T. We give the NNN negativity 
A/1,1 in Fig. 12. In the region of J2 < 0.282, there is no 
negativity Ai.i at any temperature. The threshold J2th 
is increased by the increasing temperature, similar to the 
four-spin case. 

IV. CONCLUSION 

In conclusion, we have studied the entanglement prop- 
erties of the (1/2,1) mixed-spin systems described by the 
Heisenberg model. In the systems only with NN exchange 
interactions, for two-spin and three-spin cases analytical 
results of the negativity have been obtained, which facil- 
itate our discussions of entanglement. The analytical ex- 
pression of threshold temperature after which the entan- 
glement vanishes are obtained for the two-site case. For 
the case of even number of particles, it has been found 
that the pairwise thermal entanglement is solely deter- 
mined by the internal energy, and thus builds an inter- 
esting relation between the microscopic quantity, entan- 
glement, and the macroscopic thermal dynamical func- 
tion, the internal energy in the mixed-spin systems. For 
the odd number of particles such as the three-site case, 
we also provide a relation between the internal energy 
and negativities. We have numerically studied the ef- 
fects of different finite temperature and magnetic fields 
on entanglement. As a conclusion, the thermal fluctua- 
tion suppresses the entanglement, and entanglement may 
change evidently at some critical points of magnetic field. 

In the systems also with NNN interactions, by applying 
the angular momentum coupling method, we obtained 
analytical results of all the eigenenergies of the four-spin 
system, based on which the negativity of the NN spins 
has been obtained. We also considered the excited-state 
entanglement. At finite temperature, from the partition 
function, the analytical results of negativity has been 



given. We have numerically studied the effects of the 
NNN interaction on the NN entanglement and NNN en- 
tanglement. It is natural to see that the NNN interaction 
suppresses the NN entanglement, and enhances the NNN 
entanglement. We found that the negativity between two 
spin ones is sensitive to J2 and displays some interest- 
ing properties. At finite temperature, the thermal fluc- 
tuation suppresses both the NN entanglement and the 
NNN entanglement. The threshold values J2th and Tth 
are studied in detail. The entanglement displays some 
peculiar properties, which are quite different from those 
of the spin-half model. These are due to inherent mixed- 
spin character of our system. It is more interesting to 
study entanglements in other mixed systems and explore 
some universal properties, which are under consideration. 



APPENDIX A: TWO-SITE HEISENBERG 
MODEL WITH A MAGNETIC FIELD 

The Hamiltonian of the two-site Heisenberg model 
with a magnetic field is written explicitly as 

Hi = Si^ (g) S2x + Sly (8) S2y + Si2 «) S2z + B{siz + S2z)- 

Following the same way as the discussions of subsection 
II. A, the density matrix of the thermal state is given by 
Eq. (O with the matrix elements 

ai =aee'^^ = ei^^B-^\ 

a, =a2e/^^ = ie^(e^-f2e-^), 

62=6ie''- = -f e^(e^-e-l), 

(AI) 

and the partition function 

^ ^ p|(3B-l) _j_ p~^|(3B + l) 

+ ei'^ cosh (^) cosh (^). (A2) 

Having obtained the analytical expressions of the matrix, 
we directly obtain the negativity after substituting the 
matrix elements to Eq. (|13|l . 



[1] M. A. Nielsen, Ph. D thesis, University of Mexico, 1998, 

quant-ph/0011036 
[2] M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 

87, 017901 (2001). 
[3] D. Gunlycke, V. M. Kendon, V. Vedral, and S. Bose, 

Phys. Rev. A 64, 042302 (2001). 



[4] X. Wang, Phys. Rev. A 64, 012313 (2001). 

[5] X. Wang, H. Fu, and A. I. Solomon, J. Phys. A; Math. 

Gen. 34, 11307(2001); X. Wang and K. M0lmer, Eur. 

Phys. J. D 18, 385(2002). 
[6] G. L. Kamta and A. F. Starace, Phys. Rev. Lett. 88, 

107901 (2002). 



11 



K. M. O'Connor and W. K. Wootters, Phys. Rev. A 63, [25] 
0520302 (2001). [26] 
D. A. Meyer and N. R. Wallach, J. Math. Phys. 43, 4273 
(2002). [27] 
T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 
032110 (2002). [28] 
A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 
416, 608 (2002). [29] 
Y. Sun, Y. G. Chen, and H. Chen, Phys. Rev. A 68, 
044301 (2003). 

L. F. Santos, Phys. Rev. A 67, 062306 (2003). 

Y. Yeo, Phys. Rev. A 66, 062312 (2002). 

D. V. Khveshchenko, Phys. Rev. B 68, 193307 (2003). 

L. Zhou, H. S. Song, Y. Q. Guo, and C. Li, Phys. Rev. [30] 

A 68, 024301 (2003). [31] 

G. K. Brennen, S. S. Bullock, Phys. Rev. A 70, 52303 [32] 

(2004). 

R. Xin, Z. Song, and C. P. Sun, q uant-ph/0411177| [33] 
G. Toth, quarit-ph/0406061 [34] 
F. Verstraete, M. Popp, and J. I. Cirac, Phys. Rev. Lett. 
92, 027901 (2004). 

F. Verstraete, M. A. Martin-Delgado, J. I. Cirac, Phys. [35] 
Rev. Lett. 92, 087201 (2004). 

J. Vidal, G. Palacios, and R. Mosseri, Phys. Rev. A 69, [36] 
022107 (2004). 

N. Lambert, C. Emary, and T. Brandes, Phys. Rev. Lett. [37] 
92, 073602 (2004). 

S. J. Gu, S. S. Deng, Y. Q. Li, and H. Q. Lin, Phys. Rev. [38] 
Lett. 93, 086402 (2004). [39] 
Y. Chen, P. Zanardi, Z. D. Wang, and F. C. Zhang, 
quant-ph/0407228 



S. B. Li and J. B. Xu, Phys. Lett. A 311, 313 (2003). 
K. Audenaert, J. Eisert, M. B. Plenio, and R. F. Werner, 
Phys. Rev. A 66, 042327 (2002). 

L. A. Wu, M. S. Sarandy, and D. A. Lidar, 
quant-ph/0407056 

G. Vidal, J. L Latorre, E. Rico, and A. Kitaev Phys. Rev. 
Lett. 90, 227902 (2003). 

S. Bose, Phys. Rev. Lett. 91, 207901 (2003); V. Subrah- 
manyam, Phys. Rev. A 69, 034304 (2004); M. Christandl, 
N. Datta, A. Ekert, and A. J. Landahl, Phys. Rev. Lett. 
92, 187902 (2004); Y. Li, T. Shi, Z. Song, and C. P. 
Sun, quant-ph/0406159 M. B. Plenio and F. L Semiao, 
quant-ph/0407034 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

J. Schliemann, Phys. Rev. A 68, 012309 (2003). 

L. Zhou, X. X. Yi, H. S. Song, and Y. Q. Guo, 

quant-ph/0310169 

Y. Q. Li and G. Q. Zhu, submitted. 

A. Peres, Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, 
P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 
(1996). 

G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 
(2002). 

I. Bose and E. Chattopadhyay, Phys. Rev. A 66, 062320 
(2002). 

S. J. Gu, H. B. Li, Y. Q. Li, and H. Q. Lin, Phys. Rev. 

A 70, 052302 (2004). 

X. Wang, Phys. Lett. A 329, 439 (2004). 

X. Wang and P. Zanardi, Phys. Lett. A 301, 1 (2002). 



